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Abstract 

The nanotechnology industry is growing rapidly, leading to concerns about the potential ecological consequences of the 
release of engineered nanomaterials (ENIVls) to the environment. One challenge of assessing the ecological risks of ENMs is 
the incredible diversity of ENMs currently available and the rapid pace at which new ENMs are being developed. High- 
throughput screening (HTS) is a popular approach to assessing ENM cytotoxicity that offers the opportunity to rapidly test 
in parallel a wide range of ENMs at multiple concentrations. However, current HTS approaches generally test one cell type at 
a time, which limits their ability to predict responses of complex microbial communities. In this study toxicity screening via a 
HTS platform was used in combination with next generation sequencing (NGS) to assess responses of bacterial communities 
from two aquatic habitats. Lake Michigan (LM) and the Chicago River (CR), to short-term exposure in their native waters to 
several commercial Ti02 nanomaterials under simulated solar irradiation. Results demonstrate that bacterial communities 
from LM and CR differed in their sensitivity to nano-TiOj, with the community from CR being more resistant. NGS analysis 
revealed that the composition of the bacterial communities from LM and CR were significantly altered by exposure to nano- 
Ti02, including decreases in overall bacterial diversity, decreases in the relative abundance of Actinomycetales, 
Sphingobacteriales, Limnohabitans, and Flavobacterium, and a significant increase in Llmnobacter. These results suggest 
that the release of nano-Ti02 to the environment has the potential to alter the composition of aquatic bacterial 
communities, which could have implications for the stability and function of aquatic ecosystems. The novel combination of 
HTS and NGS described in this study represents a major advance over current methods for assessing ENM ecotoxicity 
because the relative toxicities of multiple ENMs to thousands of naturally occurring bacterial species can be assessed 
simultaneously under environmentally relevant conditions. 
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Introduction 

Engineered nanomaterials (ENMs) have at least one dimen- 
sion less than 100 nm, which leads to unique chemical and 
physical properties, including high reactivity and conductivity. 
ENMs have a wide range of commercial applications, from 
pharmaceuticals and cosmetics to tools and electronics [1], and 
the nanotechnology industry has grown rapidly. For example, 
the manufacture of titanium dioxide (TiOj) nanoparticles is a 
significant niche industry and demand for these materials is 
strong. TiO, nanoparticle (nano-Ti02) production was 50,400 
tons in 2010, representing 0.7% of the overall Ti02 market. By 
2015 production of nano-Ti02 is projected to increase to 
201,500 tons [2]. Release of ENMs into the environment will be 
an inevitable consequence of high production and use [3] , with 
models predicting the entry of ENMs into aquatic ecosystems 



[4]. There is growing concern that ENMs will have adverse 
effects on aquatic ecosystems, as previous studies have 
demonstrated ENM toxicity to aquatic organisms including fish, 
invertebrates and microorganisms (for reviews see [5,6]). The 
toxicity of nano-Ti02, for example, is caused primarily by its 
photo-induced production of reactive oxygen species (ROS), 
which are powerful oxidizing agents that can damage a variety 
of cell components [7,8]. 

One of the challenges to assessing ecological risks of ENMs is 
the incredible diversity of ENMs currently available and the 
rapid pace at which new ENMs are being developed [9,10]. A 
recent survey identified 1,317 nanotechnology-based consumer 
products produced by 587 companies in 30 countries [11]. 
Assessing the relative toxicities of these materials will be 
challenging but is essential to development of ENMs with 
lower environmental impacts [12]. High-throughput screening 
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(HTS) based on robotic liquid handling and fluorescence or 
absorbance assays has been recommended as a useful tool to 
assess ENM cytotoxicity [12-14]. HTS can be used to rapidly 
test a wide range of ENMs at multiple concentrations in 
parallel, and HTS studies have provided valuable information 
on the toxicity of a variety of ENMs [15-18]. However, current 
HTS approaches for assessing ENM toxicity generally test one 
cell type at a time, with bacteria commonly used as model 
organisms [15,16,18-21]. While single-species studies can assess 
comparative toxicities of different ENMs, these studies are 
limited in their ability to predict responses of multi-species 
bacterial communities. We recently used HTS to demonstrate 
that four bacterial species common in aquatic habitats varied 
dramatically in their responses to short-term exposure to nano- 
Ti02 [22], suggesting that bacterial community responses to 
ENM exposure are complex. Therefore, the goal of the present 
study is to use HTS to assess acute responses of aquatic 
bacterial communities to ENM exposure, using nano-Ti02 as a 
model ENM. Native bacterial communities were collected from 
two aquatic habitats that differed in physicochemical character- 
istics and bacterial community composition. Lake Michigan and 
the Chicago River. Bacterial communities were exposed in their 
native waters to four different commercial nano-Ti02 materials 
in a HTS format and bacterial viability was assessed using a 
standard fluorescence assay. In addition, the effects of nano- 
Ti02 on bacterial community composition were investigated 
using next generation sequencing of bacterial 16S rRNA genes. 

Materials and Methods 

Nano-Ti02 materials 

Four nano-Ti02 materials were used: Aeroxide P25 (Evonik, 
Germany), Pigment White 6 (PW6, J.T. Baker, 4162-01), anatase 
nanopowder (ANP, Sigma, 637254), and rutile nanopowder (RNP, 
Sigma, 637262). We previously reported detailed characterizations 
of these materials [15]. Stock nano-Ti02 solutions were prepared 
at 2 g L 'in miUi-Q^ water (MUlipore, 18 Mil- cm) and sonicated 
in an ultrasonic bath at 42 kHz for 30 min. Stock solutions were 
further diluted in milli-Q_ water to the concentrations required for 
the toxicity assays. 

Sample collection 

Water samples were collected from Lake Michigan (LM) and 
the Chicago River (CR). LM water was collected in Chicago, IL 
(42.01235°N/87.66092°W) on September 27, 2012. No permit is 
required for collecting water from LM. CR water was collected in 

Morton Grove, IL (42.05904°N/087.77234°W) on October 4, 
2012. A permit for collecting water from CR was issued by the 
Forest Preserve District of Cook County. Water was collected from 
both locations in sterile Pyrex botdes, placed on ice for transport to 
the laboratory and analyzed that same day. Water samples were 
filtered using Grade 413 filter paper (VWR 28310-128). 
Characterization of these waters was reported previously [22]. 

Toxicity screening via high throughput screening 
platform 

For the HTS procedure bacterial viability was assessed using the 
Live/Dead BacLight kit (Molecular Probes). This assay allows 
determination of the relative abundance of viable bacterial cells 
based on the relative signal intensity of two fluorescent dyes: 
SYT09, which is membrane permeable, and propidium iodide 
(PI), which is not membrane permeable and quenches SYT09. A 
decrease in the ratio of fluorescent signals produced by SYT09 
(green) and PI (red) indicates a decrease in the number of viable 



bacterial cells, with viable cells defined as those with an intact ceU 
membrane. The BacLight kit has been used previously in HTS 
assays by our group and others [15,16,18,20]. 

In a preliminary experiment LM and CR water were added 
directly to microwell plates and bacterial viability was assessed 
with the Ba< Light kit. For both waters tlu; SYT09 signal was not 
above background (the signal produced by 0.2 Jt™ filtered water) 
due to the low concentration of bacterial ceUs in these natural 
waters (approximately 200 ceUs ml ' based on plate counts on 
nutrient agar). Therefore, bacterial communities in LM and CR 
water were concentrated by a fa[:tor of 200 via centrifugation at 
3,452 xg for 20 min followed by resuspension of cells in a smaller 
volume of the original source water, resulting in approximately 10* 
ceUs ml ' (based on plate counts on nutrient agar) and a SYT09 
signal that was 4 to 6 times higher than the background. This 
signal to background ratio provided an adequate dynamic range 
for assessing toxicity. 

The HTS procedure described previously [15] was used with 
some modifications. Briefly, 25 jtL of concentrated LM or CR 
bacterial communities and 25 |a,L nano-Ti02 suspensions (final 
concentrations of 0.5, 2.5, 5, 10, 25, 50 and 100 mg L~') were 
added to wells of 384-well clear-bottom microplates using a 
robotic liquid handler (Biomek FX, Beckman Coulter) with 4 
replicate weUs per treatment. Plates were illuminated with a 
xenon arc lamp (Model 6271, Newport) and incubated at room 
temperature (approximately 22°C) with shaking at 300 rpm on 
an orbital shaker. To avoid evaporation each plate was covered 
with an ultraclear film (Axygen, UC-500) that had a negligible 
effect on light intensity and spectrum distribution [15]. Nano- 
Ti02 was also tested without illumination by covering a section 
of each plate with aluminum foil. After 1 hour incubation, 
25 (iL BacLight solution was robotically added to each weU and 
plates were incubated for 15 min in the dark. Microplate reader 
(Synergy 4, BioTek) measured green (ex[:itation 485 nm and 
emission 530 nm) and red (excitation 485 nm and emission 
630 nm) fluorescence in each weU. The relative abundance of 
viable bacterial cells within each weU was expressed as a ratio of 
green to red fluorescence signals. 

Analysis of bacterial community composition 

Next generation sequencing (NGS) of 16S rRNA genes was 
used to profile bacterial communities within selected weUs 
immediately following toxicity assays. Specifically, four replicate 
wells of no treatment controls (which were not run through the 
HTS assay), and the illuminated 0 mg L ' nano-Ti02 controls, 
100 mg L~' PW6 treatments and 100 mg L~' P25 treatments for 
both LM and CR were analyzed (total of 32 samples). All material 
(75 nl) from each of these weUs was coUected and DNA was 
extracted with the WaterMaster DNA purification kit (Epicentre, 
Madison, WI). Suitability of extracted DNA for amplification was 
confirmed by PCR using bacterial primers IIF [23] and 519R 
[24] foUowed by gel electrophoresis. 

DNA from each sample was sent to Research and Testing 
Laboratory (Lubbock, TX) for tag pyrosequencing. PCR 
ampliflcation was performed using primers 28F (5'-GAGTTT- 
GATCNTGGCTCAG-3') [25] and 519R (5'-AC- 
CGCGGCTGCTGGCAC-3') [24]. Sequencing reactions uti- 
lized Roche 454 FLX (Roche, Indianapolis, IN) with Titanium 
reagents. Sequences were processed using MOTHUR (v. 1.20.8) 
[26]. Barcodes and primers were removed and low-quality 
sequences (those with an average quality score <30, any 
ambiguous characters or homopolymers >8 bases) were 
removed. Remaining high quality sequences were aligned using 
the SILVA-compatible alignment database available within 
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MOTHUR. Aligned sc-quc-iux-s \vi;rc trimmed to a uniform 
length of 250 base pairs and chimeric sequences were removed 
using UCHIME [27]. Sequences were classified by comparison 
to the MOTHUR-formatted version of the RDP training set 
(v.9) and any unknown (i.e., not identified as bacterial), 
chloroplast, mitochondrial, archaeal and eukaryotic sequences 
were removed. After pretreatment steps the LM data set 
included a total of 31,812 sequences for an average of 1,988 
sequences per sample (+/— 460) and the CR data set included a 
total of 31,253 sequences for an average of 1,953 sequences per 
sample (+/ — 253). Sequences were clustc-rc-d into operational 
taxonomic units (OTUs) based on 97% sequence identity using 
the average neighbor algorithm. For comparisons of microbial 
community composition the number of sequences within each 
OTU was used as an indicator of the abundance of that OTU 
within a sample. As with all PCR-based methods, this approach 
is subject to potential biases inherent in PGR amplification. The 
OTU abundance data were imported into Primer V.6 (Primer-E 
Ltd., Plymouth, United Kingdom) and the data for each sample 
were normalized based on the total number of sequences 
obtained for that sample. A similarity matrix was created based 
on the Bray-Curtis index, which assigns every pair of samples a 
score between 0 and 1, where 0 represents no OTUs in 
common and 1 indicates two communities containing the same 
OTUs at the same abundance [28]. An ordination of the 
samples based on the Bray-Curtis similarity matrix was 
produced via nonmetric multidimensional scaling (nMDS) [29] 
run within Primer. The analysis of similarity (ANOSIM) routine 
in Primer was used to determine if there were statistically 
significant differences in composition between the experimental 
groups [30]. Diversity of bacterial communities was quantified 
based on Shannon index [31]. 

Statistical Analyses 

One-way ANOVA was used to test the effects of nano-Ti02 
treatments on the relative abundances of viable cells, bacterial 
diversity, and relative abundances of specific bacterial phyla. 
Relative abundance data were arcsine square root transformed 
prior to ANOVA. In cases of significant treatment effects (p<0.05) 
Tukey's post hoc test was used for pairwise comparisons. 
ANOVAs and Tukey's tests were performed with Systat version 
13 (Systat Software, Inc.). 

Data sharing 

All of the sequence data analyzed in this paper can be 
downloaded from the National Center for Biotechnology Infor- 
mation (NCBI) Sequence Read Archive (SRA) with accession 
number SRP038125. 

Results 

Toxicity of nano-Ti02 Materials 

The toxicity of nano-Ti02 to bacterial communities from two 
natural surface waters, LM and CR, was assessed by HTS. These 
two waters had some differences in physicochemical characteris- 
tics: LM water had a dissolved organic carbon (DOC) concentra- 
tion of 1.56 mg L ', ionic strength of 4.75 mM and pH of 7.35, 
and CR water had a DOC concentration of 3.17 mg L ionic 
strength of 12.16 mM and a pH of 7.63 [22]. Short term (1 h) 
exposure to several forms of nano-Ti02 with simulated solar 
illumination resulted in significant decreases in the relative 
abundance of viable bacterial cells within both LM and CR as 
indicated by decreases in green to red fluorescence ratios in the 
BacLight assay. For example, exposure of LM bacteria to RNP, 



ANP, PW6 and P25 with illumination had significant effects on the 
relative abundance of viable cells at concentrations of 100, 50, 25 
and 2.5 mg L ' respectively (Figure 1). P25 showed the highi;st 
toxicity to LM bacteria, with concentrations at or above 2.5 mg 
L ' resulting in significant decreases in bacterial viability as 
compared to the 0 mg L ' control. The order of toxicity for LM 
was P25>PW6>ANP>RNP. Witiiout illumination none of die 
nano-TiOj materials had significant effects on LM bacterial 
viability (Figure 1). 

CR bacteria were less sensiti\ (: to iian<)-Ti02 than LM bacteria 
(Figure 2). For example, RNP with illumination had no significant 
effect on CR bacterial viability. ANP, PW6 and P25 with 
illumination did have significant effects on CR bacterial viability, 
although for ANP and PW6 only the highest concentration 
(100 mg L ') resulted in significantiy lower viabilities than the 
0 mg L~' control, while P25 concentrations at or above 25 mg 
L ' resulted in significant decreases in viability. The order of 
toxicity for CR was P25>PW6 = ANP>RNP. None of the nano- 
Ti02 materials tested had significant effects on CR bacterial 
viability in the absence of illumination. 

Effects of nano-Ti02 on Bacterial Community 
Composition 

Exposure of the LM bacterial community to 100 mg L ' P25 
resulted in a statistically significant shift in community 
composition as indicated by the separation of the P25 and 
0 mg L~' control samples in the nMDS ordination (Figure 3A) 
and the significant difference between these samples indicated 
by ANOSIM (Table 1). The difference in community compo- 
sition between the P25 treatment and the 0 mg L~' control was 
quite pronounced, as the average Bray-Curtis similarity between 
these treatments was only 0.14. P25 also resulted in a significant 
decrease in the di\'ersity of the LM bacterial community, as 
indicated by the Shannon index (Figure 4). Exposure of the LM 
bacterial community to 100 mg L~' PW6 also resulted in a 
statistically significant shift in composition as indicated by the 
separation of the PW6 and 0 mg L ' control samples in the 
nMDS ordination (Figure 3A) and the significant difference 
between these samples indicated by ANOSIM (Table 1). 
However, the effect of PW6 was weaker than the effect of 
P25, with the PW6 treatment and the 0 mg L ' control having 
an average Bray-Curtis similarity of 0.38. PW6 also did not 
have a significant effect on the diversity of the LM bacterial 
community (Figure 4). These results indicate that P25 exerted a 
stronger effect on the LM bacterial community than PW6, 
which is consistent with the viability results (Figure 1). Finally, 
the HTS procedure itself, i.e. exposure to simulated solar 
illumination for 1 h in microwell plates, did not have a 
significant effect on the composition of the LM bacterial 
community as indicated by the lack of separation between the 
0 mg L ' and no treatment controls in the nMDS ordination 
(Figure 3A), the lack of significant difference between these 
samples indicated by ANOSIM (Table 1), and the lack of 
significant difference in diversity between these samples 
(Figure 4). These results support the conclusion that the LM 
bacterial community changes observed in the P25 and PW6 
treatments were caused by exposure to these materials and not 
by the conditions of the HTS assay. 

NGS also revealed some significant shifts in the relative 
abundance of sequences from several of the numerically dominant 
bacterial taxa within the LM bacterial community resulting from 
acute exposure to nano-Ti02 (Figure 5; Table 2). P25 exposure 
resulted in significant decreases in sequences from the phyla 
Actinobacteria and Bacteroidetes as well as a significant decrease 
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Figure 1 . Ratio of live/dead bacteria in Lake Michigan water after 1 Kiour exposure in the dari< (D) or under simulated sunlight (L) to 
one of four nano-Ti02 materials: rutile nanopowder (RNP), anatase nanopowder (ANP), Pigment White 6 (PW) or P25. Values are 
reported as mean of 4 replicates ± standard deviation. Nt: no treatment control prior to exposure. Data points marked by asterisks are significantly 
different (p<0.05) from the 0 mg L"' nano-Ti02 control. 
doi:1 0.1 371/journal.pone.01 06280.g001 



in sequences from an unclassified bacterial phylum, with the 
Actinobacteria and Bacteroidetes sequences being almost com- 
pletely eliminated from the LM community. Within the phylum 
Actinobacteria 93% of the sequences were identified as belonging 
to the order Actinomycetales, and this order decreased signifi- 
candy in the P25 treatment. Within the phylum Bacteroidetes 
78% of the sequences were identified as belonging to the order 
Sphingobacteriales, which also decreased significantly in the P25 
treatment. As the relative abundance of these groups decreased 
with P25 treatment, the relative abundance of sequences from the 
classes Betaproteobacteria and Gammaproteobacteria increased 
within the LM communit}', although due to high variability these 
increases were not statistically significant. Within the Betaproteo- 
bacteria 57% of the sequences were identified as belonging to the 
genus Limnohacter and 36% were from the genus Limnohabitans. 
Limnobacter did show a statistically significant increase in the P25 
treatment, whereas Limnohabitans showed a significant decrease. 



Exposure of the LM bacterial community to PW6 also resulted in 
some significant changes in the composition of the community, 
including decreases in Armatimonadetes and the unclassified 
bacterial phylum, although both decreases were less pronounced 
tiian for P25 (Table 2). 

The CR bacterial community also responded to P25 treatment 
with a statistically significant shift in community composition, as 
indicated by the separation of the P25 and 0 mg L control 
samples in the nMDS ordination (Figure 3B), the significant 
difference between these samples indicated by ANOSIM 
(Table 1) and a significant decrease in diversity for the P25 
samples as compared to the 0 mg L~' control samples (Figure 4). 
The difference in composition between the P25 and 0 mg L 
tr(;atments was dramatic, with an average Bray-Curtis similarity 
of 0.12. The CR bacterial community also responded to PW6, 
with a significant difference in composition between the PW6 
treatment and the 0 mg L~' control (Table 1). However, the 
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Figure 2. Ratio of live/dead bacteria in Chicago River water after 1 Kiour exposure in the dark (D) or under simulated sunlight (L) to 
one of four nano-Ti02 materials: rutile nanopowder (RNP), anatase nanopowder (ANP), Pigment White 6 (PW) or P25. Values are 
reported as mean of 4 replicates ± standard deviation. Nt: no treatment control prior to exposure. Data points marked by asterisks are significantly 
different (p<0.05) from the 0 mg L"' nano-Ti02 control. 
doi:1 0.1 371/journal.pone.01 06280.g002 



response of the CR bacterial community to PW6 was weaker 
than the response to P25, with an average Bray-Curtis similarity 
between PW6 and 0 mg L ' treatments of 0.55 and no 
significant change in diversity (Figure 4). In fact the shift in 
CR bacterial community composition caused by PW6 exposure 
was not evident in the nMDS plot due to the much greater shift 
caused by P25 (Figure 3B); but when the CR data were re- 
plotted without the P25 samples, separation between the PW6 
and 0 mg treatments could be distinguished (Figure 3C). 
The weaker response of the CR barterial community to PW6 as 
compared to P25 is consistent with the viabihty results (Figure 2). 
Finally, the HTS procedure itself did not have a significant effect 
on the composition of the CR bacterial community as indicated 
by the lack of separation between the 0 mg L~' and no 
treatment controls in the nMDS ordination (Figure 3B and C), 
the lack of significant chflference between these samples indicated 
by ANOSIM (Table 1), and the lack of significant difference in 



diversity between these samples (Figure 4). This result again 
supports the conclusion that the changes in CR bacterial 
community composition observed in the P25 and PW6 
treatments were caused by exposure to nano-Ti02. 

NGS data demonstrated that P25 exposure resulted in some 
significant taxonomic changes in the CR bacterial community as 
compared to the 0 mg L~' control, including a significant 
increase in the relative abundance of sequences from the class 
Betaproteobacteria and significant decreases in the relative 
abundance of sequences from the phyla Actinobactcria and 
Bacteroidetes (Figure 6; Table 3). Within the Betaproteobacteria 
3 1 % of the sequences were identified as belonging to the genus 
Limnobacter and 46% were identified as belonging to the genus 
Limnohabitans , and the P25 treatments showed a significant 
increase in Limnobacter and a significant decrease in Limnoha- 
bitans as compared to the 0 mg L ' control. Within the 
Actinobacteria 95% of the sequences were identified as 
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Figure 3. nMDS analysis of composition of bacterial communities from Lake Michigan water (A) and Chicago River water (B) prior 
to treatment (no treatment) and after 1 hour exposure to 0 mg nano Ti02r 100 mg Pigment White 6 (PW6) or 100 mg 
P25 under simulated sunlight. (C) nMDS for Chicago River water with 100 mg P25 treatment excluded. Bacterial community composition was 
analyzed based on tag pyrosequencing of 16S rRNA genes. 
doi:1 0.1 371/journal.pone.01 06280.g003 



belonging to the order Actinomycetales, which showed a 

significant decrease in the P25 treatment. Finally, within the 
Bacteroidetes 62% of the sequences were identified as belonging 
to the order Sphingobacteriales and 35% to the genus 
Flavobacterium, and both of these taxa decreased significantiy 



in the P25 treatment. In contrast to P25, PW6 treatment did 

not result in dramatic changes in taxonomic composition of CR 
bacterial community, with the only significant changes being 
decreases in the order Actinomycetales and the genus Flavo- 
bacterium. 



PLOS ONE I www.plosone.org 



6 



August 2014 I Volume 9 | Issue 8 | e106280 



Effects of Nano-Ti02 on Aquatic Bacterial Communities 



Table 1. ANOSIM for bacterial community composition. 
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doi:1 0.1 371 /journal.pone.Ol 06280.t001 



Discussion 

The results of this study indicate that short-term exposure of 
bacterial communities from two aquatic habitats, Lake Michigan 
(LM) and the Chicago River (CR), to several chflFerent forms of 
nano-TiOj in their native waters under simulated solar illumina- 
tion resulted in significant decreases in the relative abundance of 
viable bacterial cells. There were some significant differences and 
some interesting similarities in the responses of the bacterial 
communities from LM and CR to nano-TiOj exposure. The 
communities from CR were generally less sensitive to nano-TiOj 
than the communities from LM, with LM communities showing 
reduced viability in response to lower doses of both PW6 and P25 
than the CR communities. This cUfference in sensitivity may have 
been due to the water matrix; we have demonstrated previously 
that when tested individually several bacterial species were more 
sensitive to nano-Ti02 in LM water than in CR water [22]. CR 



water has a higher ionic strength because of the higher level of 
dissolved salts (the concentration of chloride ions is 10 times 
greater in CR than in LM) and more importantly CR water has a 
higher concentration of DOM than LM water (the DOC of CR is 
twice that of LM), both of which can affect the behavior of nano- 
TiOj materials. For example, higher ionic strength can lead to 
increased aggregation of nano-TiOz [32], which may decrease 
production of reactive oxygen species (ROS). Our previous work 
has shown that DOM can also reduce the cytotoxicity of nano- 
Ti02 via three possible mechanisms: creating a physical barrier 
between bacteria and Ti02, absorbing UV, and ROS quenching 
[15]. These results suggest that the chemical composition of CR 
water may have contributed to lower nano-Ti02 toxicity. 
However, it is also possible that the differences between the CR 
and LM bacterial community responses to nano-Ti02 were 
related to differences in the properties of the communities 



□ Lake Michigan □ Chicago River 




No treatment OmgL-1 



PW6 



P25 



Figure 4. Diversity of bacterial communities from Lake Michigan and Chicago River prior to treatment (no treatment) and after 1 
hour exposure to 0 mg nano Ti02, 100 mg Pigment White 6 (PW6) or 100 mg P25 under simulated sunlight. Bacterial 
community composition was analyzed based on tag pyrosequencing of 16S rRNA genes. ANOVA indicated significant treatment effects for Lake 
Michigan (p<0.001) and Chicago River (p<0.05) bacterial communities. Data points marked by different letters indicate significant differences 
between nano-Ti02 treatments within Lake Michigan or Chicago River communities based on Tukey's post hoc test. 
doi:1 0.1 371/journal.pone.01 06280.g004 
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Figure 5. Mean relative abundances (in %; n = 4) of numerically dominant bacterial phyla and classes in Lake Michigan water prior 
to treatment (no treatment) and after 1 hour exposure to 0 mg nano-TiOj, 100 mg Pigment White 6 (PW6) or 100 mg 
P25 under simulated sunlight. 

doi:1 0.1 371/journal.pone.01 06280.g005 



themselves. For example, the LM community initially had a much 
higher relative abundance of an unclassified bacterial phylum, and 
this phylum proved to be especially sensitive to P25. In addition, 
previous exposure to pollutants may have played a role in CR and 
LM bacterial community responses to nano-Ti02. CR is a highly 
urbanized river with 30 combined sewer overflows (CSOs) that 



release urban runoff and untreated wastewater directly to the river 
during high rainfall events, while there are no CSO or wastewater 
effluent discharges to LM in the state of Illinois [33]. Concentra- 
tions of PCBs and chlorinated pesticides exceeding EPA standards 
have been detected in CR, with CSOs and urban runoff identified 
as likely sources, but these pollutants have not been detected in the 
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Figure 6. Mean relative abundances (in %; n = 4) of numerically dominant bacterial phyla and classes in Chicago River water prior to 

treatment (no treatment) and after 1 hour exposure to 0 mg nano-Ti02, 100 mg Pigment White 6 (PW6) or 100 mg P25 
under simulated sunlight. 

doi:10.1371/joumal.pone.0106280.g006 

water column at our LM sampling site [34] . The concentrations of suggests that ENM concentrations are likely to be higher within 

nano-Ti02 and other ENMs at these sites are unknown, as these CR than LM due to the difFerences in wastewater inputs between 

emerging pollutants are not currently monitored at these sites; but these two sites. Therefore, the high concentration of anthropo- 

wastewater has been identified as a source of ENMs [35,36], which genie pollutants in CR, possibly including nano-Ti02 and other 
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ENMs, may have selected for a bacterial community that was 
dominated by more stress tolerant taxa, resulting in a community 
that was more resistant to nano-Ti02 exposure in our study. This 
hypothesis is supported by our previous work demonstrating that 
bacterial taxa can vary dramatically in their sensitivity to acute 
nano-TiOj exposure, with some taxa showing high levels of 
resistance [22]. 

Despite the observed differences in overall sensitivity to nano- 
TiOj, when the LM and CR bacterial communities were exposed 
to a high concentration of the most toxic form of nano-Ti02 (P25) 
both showed significant shifts in the taxonomic composition of the 
communities and significant decreases in bacterial diversity as 
indicated by NGS analysis. It should be noted that the NGS assay 
is based on amplification of DNA by PGR, and it is possible to 
amplify DNA from cells that are no longer viable (e.g. that have a 
ruptured cell membrane) as long as the segment of DNA targeted 
by the PGR assay is still intact. Therefore, lack of detection of a 
cell by the NGS assay requires that the cell's DNA has been 
degraded to a point that it is not amplifiable by PGR. Previous 
studies have demonstrated the ability of nano-Ti02 to rapidly 
damage cellular DNA through the production of ROS [37]. 
Therefore, a decrease in the abundance of specific bacteria in our 
study implies not only a loss of cell viability but a degradation of 
cellular DNA. 

The specific changes in bacterial community composition were 
remarkably similar for the LM and GR communities. Specifically, 
in both the LM and GR communities P25 exposure caused 
significant decreases in the relati\'e abundances of Actinomyceta- 
les, Sphingobacteriales and Limnohahilans , and a significant 
increase in the relative abundance of Limnohacier. The CR 
community also showed a significant decrease in Flavobacterium, 
whereas Flavobacterium sequences were only occasionally detected 
within the LM community (less than 0. 1 % of the total sequences). 
Actinomycetales, Sphingobacteriales, Flavohaclerium and Limno- 
hahitans represent a phylogeneticaUy diverse collection of bacteria, 
but they have all been shown by previous studies to be abundant 
and broadly distributed in freshwater ecosystems [38-41], so their 
predominance in the LM and GR communities agrees with the 
current understanding of the ecology of bacterioplankton 
communities. The bacterial taxa Actinomycetales, Sphingobacter- 
iales, Flavobacterium and Limnohabitans all include heterotrophic 
organisms that utilize a variety of organic compounds. Bacteria 
from the genus Limnohabitans utilize algal-derived organic 
material as a substrate for growth [42]. Bacteria from the genus 
Flavobacterium are associated with the metabolism of carbohy- 
drates and polysaccharides [43] , induding cellulose [44,45] . Some 
aquatic Actinomycetes have been shown to produce cellulase and 
chitinase [46] and the Sphingobacteriales axe known for their 
ability to metabolize complex organic compounds, including 
herbicides and antimicrobial compounds [47]. Therefore, the 
observed decreases in all of these taxa with acute exposure to 
nano-Ti02 suggest that nano-Ti02 may have the potential to 
disrupt carbon cycling processes in aquatic ecosystems. 

Because the photo-toxicity of nano-Ti02 is based on the 
production of highly reactive ROS that can damage multiple cell 
structures [7,8], it is not surprising that nano-Ti02 had a negative 
effect on a wide range of bacteria taxa. What is noteworthy is that 
Limnobacler was the only bacterial taxa identified in our study that 
showed a statistically significant increase in relative abundance 
with nano-Ti02 exposure, and that this increase was observed for 
both the GR and LM communities, even though they initially had 
very different abundances of this organism. Limnobacter started 
out as one of the dominant components of the LM community, 
accounting for 27% of the total sequences, and it increased to 75% 



of total sequences after one hour of exposure to nano-Ti02. In 
contrast Limnobacter was a very minor component of the initial 
GR community (1% of total sequences) whose relative abundance 
increased dramatically (65% of the total sequences) after one hour 
of exposure to nano-Ti02. NGS data provide insight into the 
abundance of bacterial taxa based on relative rather than absolute 
gene abundance, so it is not possible to determine from these data 
if the increases in relative abundance of Limnobacter represents 
increased growth of this organism in the presence of nano-Ti02 
(i.e. the actual abundance ai Limnohacier is increasing) or simply 
higher resistan[:e to the toxic effects of nano-Ti02 (i.e. the relative 
abundance of Limnobacter is increasing because other more 
sensitive taxa are being killed). Given the short exposure time and 
the dramatic decreases in overall numbers of viable cells that we 
observed, higher resistance without additional growth seems like 
the most likely explanation for the increased relative abundance of 
Limnohacier . However, previous work by our group using the 
same HTS platform demonstrated that several bacterial species 
showed enhanced growth when exposed to nano-Ti02 in water 
from CR [22], so it is possible that the increase in relative 
abundance of Limnohacier could represent higher resistance 
coupled with enhanced growth. Whichever scenario is correct, 
our data demonstrate that Limnohacier was more resistant to the 
cytotoxic effects of nano-Ti02 than the other predominant 
bacterial groups in the LM and GR communities. The genus 
Limnobacter was first identified as a novel group of Betaproteo- 
bacteria in 2001 in a study that enriched for thiosulfate oxidizing 
bacteria from the sediment of a freshwater lake [48] and it has 
since been detected in the Baltic Sea and shown to degrade phenol 
[49]. Organisms within this genus are Gram-negative, aerobic 
heterotrophs that can grow chemolithoheterotrophically using 
thiosulfate as an energy source [48] . It is unclear at this stage why 
Limnohacier would be especially resistant to nano-Ti02, but this 
question warrants further study. The results of this study do 
demonstrate that bacterial taxa can vary significantly in their 
responses to nano-Ti02 exposure. This result is significant as it 
suggests that nano-Ti02 in the environment has the potential to 
act as a selective agent, decreasing the diversity and altering the 
composition of aquatic bacterial communities, which could have 
implications for the stability and function of aquatic ecosystems. 

Our results also demonstrate that different types of nano-Ti02 
vary in their toxicity to complex, multi-species aquatic bacterial 
communities. Specifically, the relative toxicities of the four tested 
nano-Ti02 materials were P25>PW6>ANP>RNP, which 
corresponds to the pattern we observed previously with E. coli 
[15] and to the relative photoactivity of these materials [15]. 
Since there was no effect of any of the nano-Ti02 materials on 
the CR or LM bacterial communities in the dark, our results also 
demonstrate that the toxicity of nano-Ti02 to aquatic bacterial 
communities is driven by its photoactivit}', which concurs with 
findings from previous studies using single microbial species 
[15,50,51]. 

LM and CR bacteria showed significant decreases in viability at 
P25 concentrations of 2.5 and 25 mg L ' respectively, which are 
higher than predicted concentrations for aquatic ecosystems [4]. 
However, these concentrations are relevant to possible environ- 
mental exposures. For example, Ti has been detected in the mg 
L~' range in domestic wastewater [35], and there are multiple 
scenarios under which domestic wastewater is directiy released to 
the environment. For example, in the United States untreated 
wastewater can be released to the environment via combined 
sewer overflows as well as sewer infrastructure failures or leaks. 
The U.S. EPA has predicted that the percentage of U.S. sewer 
pipes that wiU be in "poor," "very poor," or "life elapsed" 



PLOS ONE I www.plosone.org 



12 



August 2014 I Volume 9 | Issue 8 | e106280 



Effects of Nano-Ti02 on Aquatic Bacterial Communities 



conditions will increase from 23% in 2000 to 45% in 2020 [52], 
suggesting that sewer system failures and leaks are likely to increase 
in frequency in the U.S. In addition, the long-term release of low 
concentrations of nano-Ti02 to aquatic ecosystems will lead to 
accumulation of nano-TiOj in sediment, and storms or other 
physical disturbance events could resuspend nano-TiOj from 
sediments resulting in temporarily high concentrations of nano- 
Ti02 in the water column. Furthermore, the current exponential 
increase in nano-Ti02 production [53] suggests that nano-Ti02 
concentrations in the environment are likely to surpass predicted 
values. Finally, it is important to note that this study assessed only 
the potential for nano-Ti02 to kill bacteria by rupturing their cell 
membrane. It is possible that lower concentrations of nano-Ti02 
could trigger sub-lethal bacterial responses, such as changes in 
metabolic activity, gene expression, or bacterial defense systems, 
which could be relevant to the potential ecological effects of nano- 
Ti02 exposure. Further work is needed to assess these potential 
sub-lethal effects. 

One limitation of the approach described in this study was the 
need to concentrate the cells from the natural waters by 
centrifugation prior to the HTS assay. Concentrating the cells 
was necessary to achieve a cell density adequate to produce a 
fluorescent signal with the BacLight kit that was significandy above 
the background. Specifically, cells were concentrated to 10* cells 
ml This cell density is several orders of magnitude lower than 
the ceU densities generally used for HTS assays with pure cultures 
[15,16,18,22], and it is not unrealistic for natural surface waters, as 
cell densities abo\'c 10* cells ml ' have been reported for a variety 
of surface waters based on plate count assays [54—56] . In cases 
where cell densities in natural waters are above this threshold, 
concentrating cells should not be necessary to run this assay. For 
the surface waters sampled in this study, concentrating cells was 
necessary, and may have also resulted in an increased concentra- 
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might be an underestimate of nano-TiOg toxicity. 

The results of this study demonstrate that high-throughput 
screening (HTS) analysis in combination with next generation 
sequencing (NGS) can be used to assess the responses of 
complex, multi-species aquatic bacterial communities to short 
term exposure to nano-TiOj in their native waters. This novel 
approach represents a major advance over current HTS methods 
for assessing ENM ecotoxicity because the relative sensitivities of 
hundreds to thousands of naturally occurring bacterial species 
can be assessed simultaneously under environmentally relevant 
conditions. The HTS approach also offers some distinct 
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because high-density microwell plates and liquid handling robots 
enable much more highly parallel experiments than would be 
possible with other approaches, permitting multiple ENMs to be 
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